Standardized hierarchical adaptive Lp regression for noise robust focal epilepsy source reconstructions.

OBJECTIVE: To investigate the ability of standardization to reduce source localization errors and measurement noise uncertainties for hierarchical Bayesian algorithms with L1- and L2-norms as priors in electroencephalography and magnetoencephalography of focal epilepsy. METHODS: Description of the standardization methodology relying on the Hierarchical Bayesian framework, referred to as the Standardized Hierarchical Adaptive Lp-norm Regularization (SHALpR). The performance was tested using real data from two focal epilepsy patients. Simulated data that resembled the available real data was constructed for further localization and noise robustness investigation. RESULTS: The proposed algorithms were compared to their non-standardized counterparts, Standardized low-resolution brain electromagnetic tomography, Standardized Shrinking LORETA-FOCUSS, and Dynamic statistical parametric maps. Based on the simulations, the standardized Hierarchical adaptive algorithm using L2-norm was noise robust for 10 dB signal-to-noise ratio (SNR), whereas the L1-norm prior worked robustly also with 5 dB SNR. The accuracy of the standardized L1-normed methodology to localize focal activity was under 1 cm for both patients. CONCLUSIONS: Numerical results of the proposed methodology display improved localization and noise robustness. The proposed methodology also outperformed the compared methods when dealing with real data. SIGNIFICANCE: The proposed standardized methodology, especially when employing the L1-norm, could serve as a valuable assessment tool in surgical decision-making.

Pressure-Poisson equation in numerical simulation of cerebral arterial circulation and its effect on the electrical conductivity of the brain.

BACKGROUND AND OBJECTIVE: This study considers dynamic modeling of the cerebral arterial circulation and reconstructing an atlas for the electrical conductivity of the brain. Electrical conductivity is a governing parameter in several electrophysiological modalities applied in neuroscience, such as electroencephalography (EEG), transcranial electrical stimulation (tES), and electrical impedance tomography (EIT). While high-resolution 7-Tesla (T) Magnetic Resonance Imaging (MRI) data allow for reconstructing the cerebral arteries with a cross-sectional diameter larger than the voxel size, electrical conductivity cannot be directly inferred from MRI data. Brain models of electrophysiology typically associate each brain tissue compartment with a constant electrical conductivity, omitting any dynamic effects of cerebral blood circulation. Incorporating those effects poses the challenge of solving a system of incompressible Navier-Stokes equations (NSEs) in a realistic multi-compartment head model. However, using a simplified circulation model is well-motivated since, on the one hand, the complete system does not always have a numerically stable solution and, on the other hand, the full set of arteries cannot be perfectly reconstructed from the MRI data, meaning that any solution will be approximative. METHODS: We postulate that circulation in the distinguishable arteries can be estimated via the pressure-Poisson equation (PPE), which is coupled with Fick's law of diffusion for microcirculation. To establish a fluid exchange model between arteries and microarteries, a boundary condition derived from the Hagen-Poisseuille model is applied. The relationship between the estimated volumetric blood concentration and the electrical conductivity of the brain tissue is approximated through Archie's law for fluid flow in porous media. RESULTS: Through the formulation of the PPE and a set of boundary conditions (BCs) based on the Hagen-Poisseuille model, we obtained an equivalent formulation of the incompressible Stokes equation (SE). Thus, allowing effective blood pressure estimation in cerebral arteries segmented from open 7T MRI data. CONCLUSIONS: As a result of this research, we developed and built a useful modeling framework that accounts for the effects of dynamic blood flow on a novel MRI-based electrical conductivity atlas. The electrical conductivity perturbation obtained in numerical experiments has an appropriate overall match with previous studies on this subject. Further research to validate these results will be necessary.

Multi-compartment head modeling in EEG: Unstructured boundary-fitted tetra meshing with subcortical structures.

This paper introduces an automated approach for generating a finite element (FE) discretization of a multi-compartment human head model for electroencephalographic (EEG) source localization. We aim to provide an adaptable FE mesh generation tool for EEG studies. Our technique relies on recursive solid angle labeling of a surface segmentation coupled with smoothing, refinement, inflation, and optimization procedures to enhance the mesh quality. In this study, we performed numerical meshing experiments with the three-layer Ary sphere and a magnetic resonance imaging (MRI)-based multi-compartment head segmentation which incorporates a comprehensive set of subcortical brain structures. These experiments are motivated, on one hand, by the sensitivity of non-invasive subcortical source localization to modeling errors and, on the other hand, by the present lack of open EEG software pipelines to discretize all these structures. Our approach was found to successfully produce an unstructured and boundary-fitted tetrahedral mesh with a sub-one-millimeter fitting error, providing the desired accuracy for the three-dimensional anatomical details, EEG lead field matrix, and source localization. The mesh generator applied in this study has been implemented in the open MATLAB-based Zeffiro Interface toolbox for forward and inverse processing in EEG and it allows for graphics processing unit acceleration.

In silicostudy of the effects of cerebral circulation on source localization using a dynamical anatomical atlas of the human head.

Objective.This study focuses on the effects of dynamical vascular modeling on source localization errors in electroencephalography (EEG). Our aim of thisin silicostudy is to (a) find out the effects of cerebral circulation on the accuracy of EEG source localization estimates, and (b) evaluate its relevance with respect to measurement noise and interpatient variation.Approach.We employ a four-dimensional (3D + T) statistical atlas of the electrical properties of the human head with a cerebral circulation model to generate virtual patients with different cerebral circulatory conditions for EEG source localization analysis. As source reconstruction techniques, we use the linearly constraint minimum variance (LCMV) beamformer, standardized low-resolution brain electromagnetic tomography (sLORETA), and the dipole scan (DS).Main results.Results indicate that arterial blood flow affects source localization at different depths and with varying significance. The average flow rate plays an important role in source localization performance, while the pulsatility effects are very small. In cases where a personalized model of the head is available, blood circulation mismodeling causes localization errors, especially in the deep structures of the brain where the main cerebral arteries are located. When interpatient variations are considered, the results show differences up to 15 mm for sLORETA and LCMV beamformer and 10 mm for DS in the brainstem and entorhinal cortices regions. In regions far from the main arteries vessels, the discrepancies are smaller than 3 mm. When measurement noise is added and interpatient differences are considered in a deep dipolar source, the results indicate that the effects of conductivity mismatch are detectable even for moderate measurement noise. The signal-to-noise ratio limit for sLORETA and LCMV beamformer is 15 dB, while the limit is under 30 dB for DS.Significance.Localization of the brain activity via EEG constitutes an ill-posed inverse problem, where any modeling uncertainty, e.g. a slight amount of noise in the data or material parameter discrepancies, can lead to a significant deviation of the estimated activity, especially in the deep structures of the brain. Proper modeling of the conductivity distribution is necessary in order to obtain an appropriate source localization. In this study, we show that the conductivity of the deep brain structures is particularly impacted by blood flow-induced changes in conductivity because large arteries and veins access the brain through that region.

Reconstructing subcortical and cortical somatosensory activity via the RAMUS inverse source analysis technique using median nerve SEP data.

This study concerns reconstructing brain activity at various depths based on non-invasive EEG (electroencephalography) scalp measurements. We aimed at demonstrating the potential of the RAMUS (randomized multiresolution scanning) technique in localizing weakly distinguishable far-field sources in combination with coinciding cortical activity. As we have shown earlier theoretically and through simulations, RAMUS is a novel mathematical method that by employing the multigrid concept, allows marginalizing noise and depth bias effects and thus enables the recovery of both cortical and subcortical brain activity. To show this capability with experimental data, we examined the 14-30 ms post-stimulus somatosensory evoked potential (SEP) responses of human median nerve stimulation in three healthy adult subjects. We aim at reconstructing the different response components by evaluating a RAMUS-based estimate for the primary current density in the nervous tissue. We present source reconstructions obtained with RAMUS and compare them with the literature knowledge of the SEP components and the outcome of the unit-noise gain beamformer (UGNB) and standardized low-resolution brain electromagnetic tomography (sLORETA). We also analyzed the effect of the iterative alternating sequential technique, the optimization technique of RAMUS, compared to the classical minimum norm estimation (MNE) technique. Matching with our previous numerical studies, the current results suggest that RAMUS could have the potential to enhance the detection of simultaneous deep and cortical components and the distinction between the evoked sulcal and gyral activity.